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1. Introduction 

Polymers are macromolcculcs which are modelized by self-avoiding or directed 
random walks. Take for instance S = (S n ) n >o a random walk on Z starting at 
and such that |<S n +i — S n \ < 1. By polymer of dimension 1+1 and size N we 
will mean a realization of the directed random walk {(n, S n )} 0<n<N , where each 
segment [(n, S n ), (n + 1, SVi+i)] stands for a constitutive unit, called monomer. 

Suppose now that a reward h is given to a configuration {(n, S n )} 0<n<N 
each time it touches the interface, i.e. each time S n = 0. One can then consider 
a distribution on polymers of size N whose density with respect to the initial 
distribution is equal, up to a renormalizing constant, to the Boltzmann factor 

exp x Card{n 6 {1, ... , N}\S n = 0}) . 

Depending on the sign of h. this distribution favorizes or penalizes polymers 
pinned to the interface, and letting N go to infinity, the model, called homoge- 
neous pinning model, undergoes a localization/delocalization transition. 

Pinning models can also be used to study the interaction between two poly- 
mers, since the difference of two random walks is still a random walk. One can 
think for example of the two complementary strands of a DNA molecule: in this 
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case, the values of n for which S n = are the sites where the two strands are 
pinned, and the derealization transition corresponds to DNA denaturation (or 
melting). One could argue that the binding strength between the two strands 
actually depends on the base pair, which is A-T or G-C. This corresponds to 
looking at a disordered model, i.e. a model in which the reward is n-dependent. 
An assumption usually made is that the reward at site n writes 

h n = h + /3u) n 

where h € ffi, P > and ui = (w„)„>o is a frozen realization of a sequence 
of independent standard gaussian random variables. The space of parameters 
is then partitioned in localized and delocalized phases, separated by a critical 
curve (3 <— > h c (/3). The presence of disorder has important consequences on the 
model. For example, one can show that there is localization for h < provided 
that disorder is strong enough (i.e. (3 large enough) . If we consider the annealed 
model (i.e. the model in which the Boltzmann factor is averaged over disorder) , 
we have the following lower bound: 

ftcC8)>-logP(n<+oo)-^- (1) 

where T\ is the first return time of S to 0. In the last few years, many rigorous 
results were given on relevance of disorder, which in particular answer the fol- 
lowing question: when is (1) an equality? For these questions, as well as classical 
results on homogeneous and disordered pinning models, we refer to [6], [7], [13] 
and references therein. 

In this paper we remove the independence assumption on u and study the 
effect of correlations on the right-hand side of (1), i.e on the annealed critical 
curve. This is partly motivated by the long-range correlations in DNA sequence, 
see [4] and [10] on this topic. We also mention [1] and [9] where the effect of 
sequence correlation is investigated, in somewhat different contexts. In [1], the 
authors study the effect of a pulling force applied to the extremity of a DNA 
strand on the number of broken base pairs (unzipping of DNA) in two correlated 
scenarii: integrablc and nonintcgrablc correlations. In [9], the authors consider 
the effect of sequence correlation on the bubble size distribution: by bubbles 
we mean broken base pairs, and if we keep in mind the analogy with pinning 
models, it corresponds to the excursions of the directed random walk between 
two visits at 0. 

The disorder sequence in our model is a finite-order moving average of an 
i.i.d sequence, which is the simplest correlated sequence one can look at, and 
the reason for this choice will be clearer further in the text. This will be defined 
in Section 2, as well as the renewal sequence r = (r n ) n >o (the contact points) 
and the polymer measures. In Section 3, we introduce classical notions for these 
models: the free energy, the phase diagram and the (quenched and annealed) 
critical curve of the model. In the proof of Theorem 3.1, a new homogeneous 
model emerges, whose hamiltonian does not only depends on the number of 
renewal points but also on their mutual distances. Section 4 is dedicated to our 
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results and their proofs. We show in Theorem 4.1 that the difference between the 
annealed critical curve in the correlated case and the annealed critical curve in 
the i.i.d. case can be expressed in terms of the Perron- Frobenius eigenvalue of an 
explicit transfer matrix, and we give a weak disorder asymptotic of this quantity 
in Theorem 4.2. Note that the appearance of Perron-Frobenius eigenvalues is 
reminiscent of results on periodic copolymers, see [3]. 

2. The model 

2.1. Contact points between the polymer and the line 

We follow the renewal theory approach of pinning. Let r be a discrete renewal 
process such that To = and r„ = Y^—n ^fci where the inter-arrival times 
(or jumps) Tk arc i.i.d. random variables taking values in N . Furthermore, 
K (n) = P(Xb = n) = where a > and L is a slowly varying function. If 

K(oo) := 1 — J2 n >i K\ n ) > (rcsp. = 0), we say that the renewal process is 
transient (resp. recurrent). If r is recurrent, we can distinguish between positive 
recurrence (a > 1 or a — 1 and L is such that X)n>i L(n)/n < +oo) and null 
recurrence (a G [0, 1) or a = 1 and L is such that~^ n>1 L(n)/n = +oo). We 
will denote by S n the indicator of the event {n G r} = TJfc>o{ r k ~ n l so that 
if in '■= supjfc > 0|rfe < TV} is the number of renewal points before N, then 

We also suppose that for all n > 1, K(n) > (which implies aperiodicity). 
This assumption seems quite restrictive, but will be necessary in Section 4.2. If 
this condition on K were not fulfilled, we would simply have to reduce the state 
space of the matrices defined in Section 4 to {n > l\K(n) > 0} q (and assume 
that K is aperiodic). 

2.2. Finite range correlations 

Let (e n ) n gz be a collection of independent standard gaussian random variables 
(independent from r), q > 1 a fixed integer, and (ao, ■ • ■ , a q ) G W. q+1 such that 
ag + . . . + a 2 q = 1. Define the disorder sequence oj = (oJ n ) n >o by the g-order 
moving average ui n = aoe n + . . . + a q e n - q - Then w is a stationary centered 
gaussian process and its covariance function p n := Cov(wo,w„) satisfies po = 1 
and n > q => p n = 0. From now, the notations P w and will be associated to 
disorder. 

We mention that our results are adaptable to non gaussian laws, provided 
they have exponential moments. 

2.3. The quenched and annealed polymer measures 

We define the (free and constraint) quenched polymer measures, which depend 
on two parameters, the averaged pinning reward h G K and the amplitude of 
disorder /? > 0: 
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N ™ k '" = ~7a~ eX P ( ir.(fiUn + h)6 n ] (l{ a =f) +S N l {si=c} ) 

aV *N,f),h," \n=l ) 



where 



Z N,0,h,u = E ^ CX P [y^W^n + h)S^\ (l{ a =f) + <5wl{a=c})^ 

is the partition function. We also define their annealed counterparts: 

d(F®V u )* Nph 



where 



1 / N 

tin, a 



'Wn + ^)<5n (l{a=f} + ^7vl{a=c}) 



^anii , a 777 7a 

3. Generalities 

Free energy, phase diagram, and critical curve 

We give some results which are well-known for i.i.d. disorder, and which can be 
generalized to ergodic disorder (see [6, Thm 4.6, p.96]). 

Proposition 3.1. For all h € R and all (3 > 0, there exists a nonnegative 
constant F(j3, h) such that for both a = c and a = f, 

F(0,h)= lim llog^ ihiU 

Af->+oo iv 

P^-aZmosi surely and in i 1 (P u ). 

Proof. If a = c, we use the Markov property as in [6, Prop 4.2, p. 91] or [7, (3.1), 
p. 12] to write 

l°E Z N+M,/3,h,u ^ 1°6 Z U,/3,h,u + log Z N,P,h,6 M u 

where is the shift operator. We then use Kingman's subadditive theorem (see 
[12]) (w is ergodic since p n ™— ^ 0, see [5, Chp 14, §.2, Thm 2]). If a = f , we 
use the same remark as in the proof of [6, Thm 4.1, p. 94]. □ 

The phase diagram K + x R is then divided into a localized phase 

£ = {(p,h)\F(p,h)>0} 

and a delocalized one 

T> = {(P,k)\F(P i k)=0}. 

For all j3, define the critical point h c {0) '■= sup{/i G M.\F((3, h) = 0}. By convex- 
ity of F (as the limit of convex functions), T> is convex so the critical curve /3 i— ¥ 
h c (/3) is concave. Moreover, it is nonincreasing and h c (0) = — log(l — A'(oo)). 
For detailed arguments, we refer to [6]. 
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3.2. Annealed free energy and annealed bound 

The first difference that occurs when dealing with correlated disorder is that 
integrating on u> the Boltzmann factor does not yield a classical homogeneous 
model (see (3) below). As a consequence, we will need an additional argument 
to define the annealed free energy. 

Lemma 3.1 (Hammersley's approximate subadditivity [8]). Let h : N M> K be 
such that for all n, m > 1, 

h(n + m) < h(n) + h(m) + A(n + m), 

with A a non decreasing sequence satisfying: 



Mr) 
f^[r{r + l) 



E 



< oo. 



Then, lim - - exists and is finite. 

n— >+oo n 

Theorem 3.1. For all h E~$L and all j3 > 0, there exists a nonnegative constant 
F mn ((5, h) such that for both a = c and a= f, 

Moreover, if h™ n (/3) := sup{/i 6 M|F ann (/3, h) = 0} then 

KW) > hrW). (2) 

Proof. First, we compute the variance (with respect to uj) of En=i Far 
every realization of r, we have: 

(N \ N N 

u n 5 n j = cpv(uji,ujj)5iSj = s n + 2 p^SiSj. (3) 
ra=l / i,j=l ri=l l<i<j<N 

Then, = E (cxp ((/i + f ) £^ =1 <5n + /3 2 E^i* £f=7 M<fc+fc) • Now, 

we want some sort of superadditivity for the annealed partition function. For a 
polymer of size N + M, observe that 

X Pj-i S i S 3 = E + Yl Pj-i S i 6 j 

l<i<j<N+M l<i<j<N N+l<i<j<N+M 

+ Yl Pj-i^iSj. 

l<i<N<j<N+M 

Conditioned on the event {N e t}, the second term has the same law as 
12i<i<j<M Pj-i&iSj- Moreover, the third term is greater than a constant C only 
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depending only p and q. We can then write 

( ( al N + M 
>E exp (ft+y) ^ i + /3 2 Pj-iWj I SnS n + 



--1 l<i<j<N+M 



and we conclude (for a = c) by using Lemma 3.1 to — logZj^'^ with A(n) = 
—C/3 2 .For a = f , we use the same remark as in [6, Thm 4.1, p. 94] (observe that 
for a fixed value of the parameters h and /3, the amount of energy necessary to 
complete an excursion is bounded). As in [6, Prop 5.1] we use Jensen's inequality 
to prove that F(/3, h) < i ;iann (/3, h), which in turn yields the annealed bound 
(2). □ 

When disorder is i.i.d, (2) becomes h c {/3) > h c (0) - /3 2 /2 := hf 1 n (/3) and the 
question of knowing whether this is an equality was studied in several papers 
and monographs (for example, [7], [13], [6] and references therein) where we 
learn that the answer depends on the values of a and (5. 

In the next subsection, we will study the effect of correlations on /i™ n in our 
model. 

4. The annealed critical curve 

4-1- The result for q — 1 and the reason why the technique used 
does not apply to q > 1 

Proposition 4.1. If q = 1 then we have 

*™ = M0)-^-iog(i + T ^K-i)) 

Proof. If q — 1, equality (3) gives: 

( ( P 2 N ^ \ \ 

Z N*f3% = E I ex P I ( h + ~Y) lN + P 1 ^ <5 " (5 ™+l ] Sn ] ' 

The energetic contribution of a jump can only take two values: h+ (2pi + 1)/3 2 /2 
if the jump has size 1 and h + (3 2 /2 otherwise. The rest of the proof is a slight 
modification of the proof of [6, Prop 1.1], except we must consider K^ q= i^ with 
K {q=1) (l) := eP^ 2 K(l) and K {q=1) (n) := K{n) if n> 1. □ 

If q > 2, the situation is more complicated because in this case we must 
consider the energetic contribution of a q-tuple of jumps instead of one of a 
single jump. For example, if q = 2, the energetic contribution of a jump of size 1 
can be /i+(l + 2pi)/3 2 /2 or ft,+ (l + 2pi + 2p 2 )/? 2 /2, depending on the value of the 
jump just before. This idea of looking at the sequence of q-tuples of consecutive 
inter-arrival times is developed in the next section. 
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4-2. An auxiliary Markov chain and the transfer matrix 

From now we assume q > 2. We will denote by t = (io, • ■ • ,tq-i) a q-tuple in 
(N*) 9 and if (i„)„>o is a sequence, then t n := (t n , . . . ,t n+q ^i). The projection 
on the first coordinate 1 1-> iO will be denoted by ttq- Let G be a function defined 
on such g-tuples by G(t) = J2t=o Pt +...+t k > and which should be interpreted 
like this: if t is the q-tuplc of the inter-arrival times of q + 1 consecutive renewal 
points on the interface, then G(t) gives the total contribution of correlations 
between disorder at theses points. 

Notice that when we compute the value of G for some g-tuple of inter-arrival 
times, any inter arrival time strictly greater than q "does not count". To put 
it more precisely, we can consider a " cemetery state" , denoted by *, and for all 
t eN* and t we define t* := il{t< 9 } + *l{t> 9 } and = ... , t*^). Then G 
can be considered as a function of t instead of t, if we adopt the following natural 
conventions: p+ = and for all t € {1, . . . , q, *}, * + 1 = t + * = *. From now we 
will use the following notations: E = {1, . . . , q,*} and K(*) = J2 n >q K( n )- 

We will now make the following remark: the sequence of g-tuples (T„)„>o is 
a Markov chain on a countable state space, and its transition probability from 
state s = (s , . . . , to state t = (t , . . . , t q _i) writes 

9-1 

Q(s,t) := {J l {fli=ti _ l} . 

fe=i 

Note that Q is irreducible because of the positiveness of the if(n)'s. We now 
define the nonnegative matrices Qp and Q^, which will play the role of transfer 
matrices, by 

n ( c 7\ - /G(i)_^M.TT -, 

°^ (M) - e l-if(oo) 11 

and 

Q " (s '* ] ~ e i-if(oo) IPK^I-J- 

v 7 i— 1 

Since is an irreducible nonnegative matrix on the finite state space E q , 
we know by the Perron-Frobenius theorem that there exists a Perron-Frobenius 
eigenvalue X(/3) and an associated right eigenvector v* = {v*(t )) with positive 
components (see [11]). 

4-3. Statement of the results 

We are now ready to state our main results. The first one expresses the annealed 
critical curve in terms of the Perron-Frobenius eigenvalue of the transfer matrix 
Q}. 
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Theorem 4.1. For all a > 0, for all f3 > 0, 

K nn (P) = Mo) - MP) 

where 

A(/3) :=^+logA(/3). 

It seems difficult to give a nice explicit expression of A(/3), since it is the 
Perron- Frobenius eigenvalue of a matrix of size (q + l) q . For q — 2, we have 
computed 

A(/?) = ^+log </>(/?) +log< V 



2 



where 



(/ 3) = i + -^-M p ^ )02 - 1) + t-^^ 2 - 1) 

1 — K(po) 1 — a(oo) 

^ 1-K(oo)\ l-K(oo)J 1 j 

V 1 - ^(oo) - K(l) y ' 

In the general case, the asymptotic behaviour of the annealed critical curve 
for weak disorder can be explicitcd: 

Proposition 4.2. Let f be the renewal process starting from with inter-arrival 
time law K := -, — — r. Then we have 

1 — K (oo) 




By a classical argument (see [7, p. 6] or [6, Rmk 1.19, p. 40]) we will restrict 
ourselves in what follows to r recurrent, i.e. if (oo) = 0, h c (0) = and f = r. 

Before going into details, we outline the proof of Theorem 4.1. First, we 
introduce in Lemma 4.1 new Markov transition kernels built from the transfer 
matrices and an eigenvector associated to A(/3). From these we give a new law for 
the sequence of q-tuples of consecutive inter-arrival times, to which we associate 
what could be called a "q-correlated" renewal process. This process is in fact 
a particular Markov renewal process (these are processes in which the return 
times are not necessarily i.i.d., but driven by a Markov chain, see [2] on this 
subject). With Lemma 4.4, we link the annealed free energy of our initial model 
to the homogeneous free energy of the new "q-correlated" renewal process. This 
will be the starting point of the proof of Theorem 4.1. Note that for positive 
recurrent renewal processes we give a shorter proof than in the general case. 
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4-4- A "q- correlated" renewal process related to the model 

For all q-tuplcs t, define u(t) = v*(t ). 

Lemma 4.1. Qp(s,t) := Q ^^ ) and Q* p (s* ,T) ■= Q ^(p)^') ^ are Markov 
transition kernels. 

Proof. For Qp, the result is a direct consequence of the relation Q*aV* = A(/3)z/* 
and of the positiveness of A(/3) and v. For Qp, we write for all s = (sq, . . . , s g _i), 

E Qp{s,t)v(t) = Y, ef G ^-' s «-^K{t)v{s u . . . , 
t *>i 

= E e^ s ''-'<-^K(t)u*(s* v . . . , s* q _ v t*) 
*>i 

= E eP 2Gi - s i'-^' t ">K{t*)v*{sl 1 ... 1 s* q _ 1 ,t*) 

f£E 

= \(j3)v*(s*) 

= WW*)- 

The result follows in the same way as for Q*p. □ 

Since Qp is a finite irreducible transition matrix (it has the same incidence 
matrix as Qp, which is irreducible), it has a unique invariant probability mea- 
sure that we denote by n* . If we define fi a measure on (N*) 9 by = 

K(to)_ Bh=ll u *(t*) then 
K(t*) ■ ■ ■ WttfZJ^ ^ >■> 

Lemma 4.2. /i is the invariant probability ofQp. 

Proof. By a direct computation of X/t M^)' fi is a probability. Now we prove 
that it is invariant. For all t £ (N*) q , we have 

Y^)Qp(s, t) 

8 

= \w-i<f°® v ®K&-x) E I s '! ' • • • f ~ 2 1 

^ V\S,to, . . .,t q -2) 

K(s)K(t ) . . . K(t q - 2 ) M*0*, to, • • • , t*_ 2 ) 



x 

s> 



^ K(s*)K(t* ) . ..K(t*_ 2 ) v*(s*,t* , ■ ■ -X q -2) 

(m-^^)K { t U) ^ e w^'Ir-'^i 
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where for the last equality we use the fact that p.* is the invariant probability 
ofQJ. □ 

Let (Tk)k>o be a process of inter-arrival times in N* defined by the ini- 
tial distribution P(T = t ,. . . , T g _ 1 = = n&=o T^kToo) anc ^ ^he re l a ti° n 
P(T fc = t q \T k -q = t ,...,T k -i = t q -i) := Qp((to, . . . ,t q -x),{tx, . . . ,t q -x,tq)) 
for all k > q. From this we define the process f by fo = and for all n > 1, 
f n = Y^f—o^k- We generalize the definitions and notations of Section 2.1 to 
the process f: P, iff, and so on. The process f is what we have called a "q- 
correlated" renewal process, since the law of one inter-arrival time depends on 
the value of the q previous ones. 

Lemma 4.3. For all h el and all j3 > 0, 

F mn (p,h) = lim llogE^+T^+^E^G^A 

fi/ie sum in the right hand side is void on the event {in = 0}/ 

Proof. On one hand, we have by integrating over disorder the partition function: 

Z^p, h = IE | cxp J i N {h +^) + f3 2 J2 

y y i<«<i<jv 

= E ^cxp ^i N {h + y) + ^^ PkSi*H-kj ^ • 

On the other hand, Y^Li G(T n ) = Y,iLi YX=i Pk^A+k- We prove the lemma 
by showing that there exists a constant C(p,q) such that 



i N N-lN-i 

E G(T n ) -EE PkSi^-k 

n—1 i—1 k—1 



<C{p,q). 



Indeed. 



N q N-1 q q 

E E Pk^i+k = E E Pk^i+k + E Pk$N&N+k 
i=l fc=l i=l fc=l fc=l 

N—lN-i q 

= E E Pk^i+k + J~] Pk&N&N+k 
i=l k=l k=l 
N-l q 

+ E E PkSiSi+k 

i=N-q+l k=N-i+l 

where the second term is bounded in absolute value by q x maXj = i... Q \pi\ and 
the third term by q< - q ^ x maXj = i... Q \pi\. □ 
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Lemma 4.4. For all h £M and all j3 > 0, 

F ann (B,h)= lim 4 logE f e (^+AW)w 

Proof. Remember that we suppose _fT(oo) = 0. By decomposing on the possible 
values of in, we have on one hand: 



N 



n=l 

J2 e? 2 ES=i G( ^Q(t , ii) . . . Q(i n _x,? n )^(i ) 



t ,...t n 

t + ...+t„_i<JV 
t + ...+t„>Ar 



e (fe+AG8))n 

E^) + E Af g)n E <9^(*0,il)...^(in-l,in)^ 9 (t0) 



^ A(/3) 

t + ...+t„_ 1 <N 
t +...+t n >N 

N 



K (t) + E e (h+Amn 



t>N 

X 



to, 

t + ...+«n-l<iV 

to + ...+*, 1 >A r 



and on the other hand, 



N 

= K (t) + E e (h+A(/3))n £ Q/3(t ,tl) ■ • ■ Q^n-l,t n )K^(t ] 



t>N n=l to,...t„ 

t +...+t n -i<W 
t + ...+i„>iV 



bmce v (?) = ) and ^* is a finite vector with positive components, there 
exists c and C two positive constants such that for all to,t n , c < v %°\ < C. We 
conclude by using this remark and Lemma 4.3. □ 



4-5. A short proof of Theorem 4-1 in the positive recurrent case 

In accordance with Lemma 4.4, we will work on the homogeneous pinning model 
of the process f. In the positive recurrent renewal-type lemma is obtained, 

which allows us to conclude. 
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Lemma 4.5. If a > 1, or if u~\ and L is such that £n>i L(n)/n < oo then 
tends almost surely and in L 1 to a positive constant. 

Proof. From Lemma 4.2, the sequence of g-tuples (Tfc, . . . , Tk+ q -i)k>o is a pos- 
itive recurrent Markov chain, with invariant probability measure fi. If the pre- 
vious conditions on a are satisfied, t:q : t — > to (the projection on the first 
coordinate) is /x-integrable. As a consequence, 

f N 1 N ~ 1 

-J7 = ^2 7r o(r fe ,...,f fe+g _i) "4 c := ^ t X ^(t 0) ...,i 9 _i) < oo. 

fc=0 fi,...,t a _i>l 

We deduce from this that ^ °4' i > by using the inequality T t ~ N 

The convergence in L 1 follows from the Dominated Convergence Theorem. □ 

From Lemma 4.4, h < -A(/3) implies that F alLn (f3, h) = 0. Suppose now that 
h = — A(/3) + e with e > 0. By Jensen's inequality, we have that 

llogE(> +A ^~) >e^^. 
We conclude that F ann (/3, h) > by using Lemma 4.5 and Lemma 4.4. 



4-6. Proof of Theorem 4-1 in the general case 

We now give a proof without any assumption on a. The starting point is Lemma 
4.4 and we will actually identify the free energy of the pinning model associated 
to f . Let's fix e > 0. We introduce the matrices 

Qp, F (s,t) =e e - Ft - l QpfaI) 

and 

Q^ F {TX) = e e - F4, ^-^Q}{rX) 
where 4>f{s*) = s* if s* 6 {1, . . . , q} and 

M*) = - ? iog-^) 

i.e. 0_f(*) verifies 

e-F.M*)^) = ^ e - Ft K(t). (4) 
We will denote by Xp the Perron- Frobenius eigenvalue of the irreducible matrix 

Lemma 4.6. There is a unique positive real denoted by F(e) such that 
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Proof. Elcmcntwise, Qp F is smooth and strictly decreasing with respect to F. 
Since Xf is a simple root of the characteristic equation of Q*^ F (see [11, Thm 
1.1]), Xf is also a smooth function of F by the Implicit Function Theorem. From 
the formula (see [11]) 

Xf = max mm — — 

u>0 j:vj>0 Vj 
Y.E1 v i = 1 

one also obtains (see [6, Appendix A. 8]) that Xf is strictly decreasing in F and 
that Xf —> as F — > oo. Since Ao = exp e > 1, the result follows. □ 

Let v* be a Perron- Frobenius right eigenvector of Q* p ,. . We define v by 

m =**(?)■ (5) 

Lemma 4.7. The matrices 

P(s,t):=Q 0iP(e) (s,t)^ (6) 

and 

p*(n**):=Q;, #(e) (^r)||l oo 

are stochastic and irreducible matrices. Furthermore, if we denote by I* the 
invariant probability measure of P* , then I defined by 

l{t):=r{T)T\ K{Sj)e — (8) 

is the invariant probability measure of P. 

Proof. The proof is left to the reader. It consists in straightforward computations 
very similar to Lemma 4.1 and Lemma 4.2. We use Lemma 4.6 to prove (7), and 
(4), (5), (7) to prove (6). □ 

Note that, like Qp and Qp, P satisfies the "consistancy " condition 
P(s,t) ^ Si = U-iVi e 1, ...,<?- l. 
We define the process by = and 

n— 1 71—1 



where (Ti F) ) „> is a Markov chain of transition P and initial distribution I. 
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Lemma 4.8. There exists two constants C > c > such that 

ce^ N P(N G t< f >) < E (exp(« w )^) < Ce^ N ¥(N G r^) 

Proof. Decomposing the partition function and using (6) we get 
E (cxp(ei N )S N ^j 

N 

= £ £ e era Q^(t 0j ti)...^(in- g -i,I„- g )if® 9 (io) 



n=l i ,...,t n -l 

t + ...+t n - g =N 

N 



= e F(e)NJ2 P(h,h)---P(tn- q -l,t n - q )l(i ) 

t + ...+tn-l=JV 

^ / £>(I ) X^(t ) e -^( e )(«Q+-+« g -i) \ 
X i(*o) 6 ) 

and, from (5), (8) and the finiteness of E q , the term in parenthesis is uniformly 
bounded by two positive constants C and c. □ 

From this we deduce: 

Lemma 4.9. For all e > 0, 

lim logE (exp(e»jv)5jv) = ^0) > 0. 

Proof. Since P(7V G t( f ') < 1, it will be sufficient to prove that 

liminfPfiV G r (F) ) > 0. 

N-Hx> 

We use an argument that has been already used in the study of Markov renewal 
processes arising in the study of periodic pinning (see [2, Chp VII. 4], [3] or 
[6, Chp 3]). We choose arbitrarily the state 11 = (1, . . . , 1) G (N*) 9 . Consider 
(@n)n>o the following sequence of stopping times: 

00 = inf{n > Q\T ( P = 11} 

B k+1 = inf{n > 9 k \T [ P = 1} 

Since (T n ) n >o is positive recurrent, these stopping times are finite almost 
surely. If we now define the process t by t® := Tq then it is clear that 

F(N G t (f) ) > ¥(N G t 9 ) 
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By the strong Markov property, r e is a (delayed) renewal process whose inter- 
arrival times are on average equal to 



m := E U (T (F) + . . . + T 6a ^) F = E u <T ir (T { n F) )1 {f>0>n} 

n=0 



^ 7 ro(i)E ]i ^l { 



n=0 



5< Jr) =t,9o>n} 



t 



< oo. 



By the Renewal Theorem, we have 



V(N € t 6 ) "4°° l/m> 

and the proof is complete. □ 

Theorem 4.1 is now a direct consequence of Lemma 4.4 combined with Lemma 
4.9. 



4-7. The weak disorder asymptotic: proof of Proposition 4-2 

We now give some lemmas which will be useful for the proof of Proposition 4.2. 
If I C E q then we will denote by Qq' 1 the matrix with entries Qq' 1 (s* ,t*) = 
Qq(s* ,t )lr£* eJ i. If P is an n by n matrix then Com(P) is the matrix of the 
cofactors of P, i.e. Com(P)(i, j) = (— dctMjj where Mij is the n — 1 by 
n — 1 matrix obtained by deleting the i-th line and the j-th column of P. 

Lemma 4.10. Qq is primitive and its invariant probability measure is K® q (s*) = 
K(sl)...K(s' q _ l ). 

Proof. For all T E E q , 

(K®«Q* )(T)= Y, K®«(s*)Q* (s*,T) 

= Y K^( s *,ti...,t* g _ 2 )K(t;_ 1 ) 

s*£E 

= Y K(s*)K(t* )...K(t* g _ 2 )K(t* g _ 1 ) 

s*£E 

= K® q {T) 
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so K® q is the invariant probability measure of Qq. Moreover, 

(Qo) q ( s ) = F (Tg = t , . . . ,T 2q _ 1 = t q _ 1 \T = Sq, . . . ,T*_! = s g _i) 
= F(T;=t* ,...,T* q _ 1 =t q _ 1 ) 
= K® q {T) 

which is positive under the assumptions of Section 2.1. Since (Qq) 9 > 0, Qq is 
primitive. □ 

Lemma 4.11. Tr (Com^ - Id)) + and for all s* G E q 

Tr (* Com(Q^Id)Q*' { ^ }N 



Tr(' Com(Q* - Id)) 



K® q (r). 



Proof. In this proof we will use the properties of the Perron- Frobenius eigenvalue 
of a primitive matrix, that one can find for example in [11]. 
We define for all s* e E q : 



p(s*) 



Tr(* Com(Q5-Id)Q;' {s * } 
Tr (* Com(Q* - Id)) 



By Lemma 4.10, we only need to prove that p is the invariant probability measure 
of Q* . 

Since Qq is stochastic, 1 is clearly a right eigenvalue of Qq with associated 
eigenvector 1L (the vector with 1 on all its components). Moreover, Qq is primi- 
tive (Lemma (4.10)) so the Pcrron-Frobenius eigenvalue exists and all we have 
to prove is that |A| < 1 for every (possibly complex) eigenvalue of Qq. Indeed, 
if v is an eigenvector associated with such an eigenvalue, and s* G E q is such 
that v(s*) = max \v(t )| then 

t*£El 

Xv(s*)= Qo(s*X)v(i*) 

so A||i>(s*) < |w(s*)|, i.e |A| < 1. This proves that 1 is the Perron- Frobenius 
eigenvalue of Qq, with associated eigenspace Ri. 

Now, from [11, Ch. 1, Corollary 2], we have that the rows of * Com(<5Q — Id) 
are all equal to the same left eigenvector (for the eigenvalue 1) of Qq, that we 
will denote by L. A first consequence is that Tr (Coiti(Qq — Id)) ^ because 
the entries of L are either all positive or all negative. Another consequence is 
that if wc define 

m = (m(s*))s* eE « = (Tr (* Com(Q* - ld)Q* Q ^ } ) ) 

then m(s*) — L(~s*) for all s* G E q . Moreover, from the relation 

(Q*q - Id)' Com(Q* - Id) = 
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we deduce that m(s*) = Tr (* Com(QQ — Id)). Since p is simply m renor- 

malized by X)s*e_E9 to (**)j it i s the invariant probability of Qg. □ 
Proof of Proposition In what follows, we will use the notations Q' and Qq 
as shortcuts for -g^-(O) and g^/ (0). First we will show that /3 H> \([3) is 
infinitely differentiable (C 2 would be enough). Let's define 4>(/3,\) = dct(Q^ — 
Aid) so that </>({3,X) is the characteristic polynomial of Q*^, and 4>(f3, A(/3)) = 
for all /3. The Perron- Frobenius eigenvalue of a nonnegative primitive matrix 
being a simple root of its characteristic equation, §f(/3, A(/3)) 7^ for all /? > 0. 
Since is infinitely differentiable, the same holds for A by the Implicit Function 
Theorem. 

Now, a straightforward computation shows that (we use that A(0) = 1) 

>> - y <°> 

|^(0) = 1+A"(0)-A'(0)'- 
All we need to show then is 



A'(0) = (9) 



9 



A"(0) = 2^p„P(T ! er). (10) 



n=l 



By derivating the relation ^>(/3, A(/3)) = we obtain 

|(0 ) l) + A'(0)g(0,l ) = 0. 

We already know that |f (0, 1)^0 and since Q' = then §|(0, 1) = 0, which 
leads to (9). 

All we have to do now is to prove (10). A Taylor expansion of det(Q^ — A(/3) Id) 
gives: 

det(Q; - AGS) Id) = det (q* -Id+(Q{( - A"(0) Id)£- + o(/? 2 )) 

= Tr ( f Com(Q5 - Id)(Q ' - A"(0) Id)) ^ + o(/3 2 ) 

where we have used the differential of the determinant: det(A + H) — det (A) + 
Tr(?Gom{A)H) + o(\\H\\). But since det(Q£ - A(/3)Id) = we have 

Tr (* Com(QS - Id)(Q ' - A"(0) Id)) = 

which yields 

= TfCom(Qg - Id)Q(Q 
1 ' Tr( t Com(Q5 - Id)) ' 
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Note that Tr(*Com(Q5 - Id)) ^ (Lemma 4.10). 

Let's now consider Qq as a function of {p n )i<n<q- Observe that 

QZ(r,r) = 2G(T)Q*(s*X) 

so Qq linearly depends on (p n )i<n<q- We have then 

Qo(f>i,...,p q ) = Qo (0, • ■ • , 0, p q ) + Q' \ Pl , . . . , p q -!, 0) 

The result of the theorem is clearly true for q = 1 (remember that we have an 
explicit expression of A in this case, see Proposition 4.1 ) so we can suppose 
that it is true for a (q — l)-order moving average and show that the result holds 
for q. The induction hypothesis then implies 

Tr(*Com(QS - Id)Qg( P i, . . . , Pq -u 0)) _ - p , c , 

Tr(*Com(Q5 - Id)) ~ ^ ^ PnFl " G Tj 



n=l 



so the only thing left to prove is that 

Tr( t Com(QS-Id)Q[ ) '(0,...,0,p g )) 



2p q V(q G r). (11) 



Tr('Com(Q* - Id)) 

Let's define I q = {s* G i? 9 s.t. p q appears in G(s*)} and notice that 

Q'^0,...,0,p q ) = 2p q Q* ' 1 ". 

We obtain from Lemma 4.11: 

TrCCom(QS-Id)Q^) = v 8 * 
Tr( t Com(Q5 - Id)) 1 j 

= P(g G t) 

which proves (11). □ 



5. Open questions 

Here are possible developments and open questions we wish to discuss. 

1. What happens to these results if we consider a more general moving av- 
erage ui n — J2k>o a k £ n~k> w hcre J2k>o a k = 1' For example, it would 
be interesting to consider power-like decaying correlations: the difficulty 
is that much of the arguments used for finite order moving averages, in 
particular the reduction to a finite space E = {l,...,g,*}, cannot be 
applied. 

2. In which cases do we have h c (/3) = /i^ nn ( / 5)? In the case of i.i.d. disorder, it 
has been shown that for a < 1/2, disorder is relevant (h c (j3) = /i* nn (/3) for 
small enough (3), while for a > 1/2, it is irrelevant. Does the same critical 
value a c = 1/2 hold in the correlated case? 



Julien Pois at /Pinning with correlated disorder 



19 



References 

[1] Allahverdyan, A.E., Gevorkian, Z.S., Hu, C.K., Wu, M.C.: Unzipping of 
dna with correlated base sequence. Phys. Rev. E 69(6), 061,908 (2004). 

[2] Asmussen, S.: Applied probability and queues, Applications of Mathemat- 
ics (New York), vol. 51, second edn. Springer- Verlag, New York (2003). 
Stochastic Modelling and Applied Probability 

[3] Caravenna, F., Giacomin, G., Zambotti, L.: A renewal theory approach to 
periodic copolymers with adsorption. Ann. Appl. Probab. 17(4), 1362-1398 
(2007) 

[4] Chen, X.Y., Bao, L.J., Mo, J.Y., Wang, Y.: Characterizing long-range cor- 
relation properties in nucleotide sequences. Chinese Chemical Letters Vol. 
14 14(5), 503-504 (2003) 

[5] Cornfeld, LP., Fomin, S.V., Sinai, Y.G.: Ergodic theory, Grundlehren der 
Mathematischen Wissenschaften [Fundamental Principles of Mathematical 
Sciences], vol. 245. Springer- Verlag, New York (1982) 

[6] Giacomin, G.: Random polymer models. Imperial College Press, London 
(2007) 

[7] Giacomin, G.: Renewal sequences, disordered potentials, and pinning phe- 
nomena (2008). URL arXiv.org: 0807. 4285 

[8] Hammcrsley, J.M.: Generalization of the fundamental theorem on sub- 
additive functions. Proc. Cambridge Philos. Soc. 58, 235-238 (1962) 

[9] Jeon, J.H., Park, P.J., Sung, W.: The effect of sequence correlation on 
bubble statistics in double-stranded dna. The Journal of Chemical Physics 
125 (2006) 

[10] Peng, C.K., Buldyrev, S.V., Goldbergcr, A.L., Havlin, S., Sciortino, F., 

Simons, M., Stanley, H.E.: Long-range correlations in nucleotide sequences. 

Nature 356, 168-170 (1992) 
[11] Seneta, E.: Non- negative matrices and Markov chains. Springer Series in 

Statistics. Springer, New York (2006) 
[12] Steele, J.M.: Kingman's subadditive ergodic theorem. Ann. Inst. H. 

Poincare Probab. Statist. 25(1), 93-98 (1989) 
[13] Toninelli, F.L.: Localization transition in disordered pinning models. In: 

Methods of Contemporary Mathematical Statistical Physics, Lecture Notes 

in Mathematics, pp. 129-176 (2009) 



